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Abstract. The nonlinear interaction between electromagnetic, electrostatic and 
gravitational waves in a Vlasov plasma is reconsidered. By using a orthonormal 
tetrad description the three-wave coupling coefficients are computed. Comparing 
q ■ with previous results, it is found that the present theory leads to algebraic expres- 

sion that are much reduced, as compared to those computed using a coordinate 
frame formalism. Furthermore, here we calculate the back-reaction on the gravita- 
bJ). tional waves, and a simple energy conservation law is deduced in the limit of a cold 

plasma. 
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On 

1. Introduction 

Much work has been devoted to special relativistic effects of plasmas [1], largely 
, stimulated by the rapid progress in high power laser technology [2] . Several studies 

of plasmas in a general relativistic context have also been made, see .e.g. [3-5,11-13], 
considering for example plasmas in strongly curved space-times close to pulsars [14], 
multi-fluid plasma effects in general relativity [15-17], dynamo effects in strong 
gravity [18], the effect of two-temperature systems on scalar perturbations [19], 
viscous heating of accretion discs due to gravitational wave dissipation [20], or 
considering gravitational wave (GW) propagation in a plasma medium, in particular 
nonlinear interactions in e.g. dusty plasmas [21,22] or MHD plasmas [23-27]. GW:s 
are currently opening up a promising new window for astronomy and astrophysics, 
e.g. astroseismology [28,29], and the interaction between GW:s and electromagnetic 
fields have been proposed as a possible means to detect GW:s [30-34]. Moreover, 
the nonlinear interaction of a curved spacetime with and electromagnetic field can 
yield a multitude of interesting astrophysical and cosmological effects [35-38]. In 
a magnetized plasma, a basic effect is the linear coupling between electromagnetic 
(EM) waves and GW:s that occurs for propagation across a static magnetic field 
[5].f This linear mechanism leads to the excitation of magnetohydrodynamic waves 
in a plasma [11-13]. Naturally linear coupling mechanisms is not sufficient if one 
is interested in the possibility of converting GW energy to frequencies different 
from that of the original source. Numerous examples of such mechanisms exists in 

f For a discussion concerning the case of cosmological magnetic fields, see Refs. [6-10]. 
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plasmas, involving e.g. frequency up-conversion due to nonlinear wave steepening 
[39], or various three- wave couplings between GWs and electromagnetic waves, e.g. 
[40,41]. Wave coupling mechanisms involving GWs are studied for several different 
reasons. In some cases, the emphasis is on the basic theory [42-45]. In other works, 
the focus is on GW detectors [46-48], on cosmology [49-51], or on astrophysical 
applications such as binary mergers [52], gamma ray bursts [53], pulsars [54] or 
supernovas [55]. 

In the present paper we will re-consider the problem of three-wave interacting 
between electrostatic (ES), electromagnetic (EM) and GW:s in a plasma, using a 
collisionfree kinetic description, i.e. the Vlasov equation [40]. In contrast to pre- 
vious authors we will use a tetrad formalism [56], rather than a coordinate frame 
formalism, since the former formalism has been shown to significantly reduce the 
algebraic complexity. In particular, the coordinate frame formalism applied on the 
present problem produced very cumbersome algebraic expressions for the coupling 
coefficients [40], and as a result it was impossible to see if the growth rates for 
parametric processes where positive in general (when damping was omitted) , which 
is related to the fulfillment of the so called Manley-Rowe relations [57]. The coeffi- 
cients derived using the present tetrad formalism is shown to agree with previous 
works in the limit of a cold plasma [40]. For the case of a finite temperature, the 
coefficients found from the present formalism are algebraically much simpler than 
previous results. However, due to the complexity of the previously computed coef- 
ficients, a comparison cannot be easily made for the general case. Finally, in the 
present paper we also include the back-reaction on the GW, in contrast to previous 
works. This allows us to discuss the energy conservation properties, and an energy 
conservation law is presented in the low temperature limit. 



2. Basic equations in the tetrad formalism 

We consider the interaction between weak gravitational waves and a collisionless 
plasma in an external magnetic field. Since we consider nonempty space the back- 
ground space-time is necessarily curved. However, if the wavelength of the gravita- 
tional waves and the interaction region is small relative to the background curvature 
we may take the background to be flat and static (the lowest order version of the 
high-frequency approximation) , and consider the perturbed energy- momentum ten- 
sor corresponding to the perturbations of the electromagnetic and material fields. 

For simplicity, the unperturbed plasma is assumed to be static, isotropic and 
homogeneous. Linearized, the Einstein field equations (EFE) take the form 

Dh ab = -2k [ST ab - \5Tn ab ] (2.1) 

provided the gauge condition h ab b = is fulfilled, which is equivalent to state that 
only tensorial perturbations are present. Here □ = [c~ 2 9 t 2 — <9 2 ] , h ab is the small 
deviation from the Minkowski background metric, i.e. g ab — n ab + h ab , n = 87rG/c 4 , 
5T ab is the part of the energy-momentum tensor containing small electromagnetic 
and material field perturbations associated with the gravitational waves and 5T 
= 5T a a . In the following it is understood that we neglect contributions of second 
order and higher in h ab . In our notations a, b, c, ... = 0,1,2,3 and i,j,k,... — 1,2,3 
and the metric has the signature ( — h H — h). 

In vacuum, a linearized gravitational wave can be transformed into the transverse 
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and traceless (TT) gauge. Then we have the following line-element and correspond- 
ing orthonormal frame basis 

ds 2 = -c 2 dt 2 + (1 + h + (0) dx 2 + (1 - h + (0) dy 2 



+2/i x (0 dxdy + dz 2 



e = c 1 d t , ei = (1 - 
e 2 = (1 + \h+)d y - \h x d.. 



\h+)d x 



\h x d, 



e 3 = d z 



(2.2) 

(2.3) 
(2.4) 



where £ = z — ct and h + ,h x <C 1. As it turns out, the gravitational waves takes 
this form also in the particular case (propagation in an isotropic plasma) that we 
are considering. The difference to the vacuum case will be that £ — z — v p ht, where 
v p h is the phase velocity of the gravitational wave. Note, however that the theory 
will be limited to the case w p h « c, due to the omittance of back-ground curvature 
effects. From now on we will refer to tetrad components rather than coordinate 
components. 

We follow the covariant approach presented in [56]] for splitting the electromag- 
netic and material fields in a 1 + 3 fashion. Suppose an observer moves with 4- 
velocity u a . This observer will measure the electric and magnetic fields E a = F a i,u b 
and B a = i e a bcF bc , respectively, where F a b is the electromagnetic field tensor and 
e a t c is the volume element on hyper-surfaces orthogonal to u a . ft is convenient to 
introduce a 3- vector notation E = (E % ) = (E 1 ^ 2 , E 3 ) etc. and V = e». From now 
on we will assume that u° — c is the only nonzero component of u a . As has been 
presented in e.g. Refs. [58, 59]] the Maxwell equations contain terms coupling the 
electromagnetic field to the gravitational radiation field. Including terms that are 
linear in h+ and h x , but omitting terms that are quadratic and of higher order, 
Maxwell's equations are written as 



V x B 

d t B 
V B 
V-E 



Mo(j +3 E ) + 
V xB-jB, 


Pc 



1 



(2.5) 

(2.6) 
(2.7) 
(2.8) 



where 



IE = -TT 



JB = 



f 

2~c 

l r 

2 
f 

2 
1 

2 



(E x - cB y )h+ + (E y + cB x )h 
(E y + cB x )h+ + (E x - cB y )h x 

(Ey + cB x )h+ - (E x - cB y )h x 

(E X - CBy)h+ + (Ey + cB x )h x 



ei 



e 2 , 



ei 



e2, 



(2.9) 



(2.10) 



are effective currents due to the GWs, see e.g. Ref. [45]], and the dot denotes 
derivatives with respect to the argument. The physical current- and charge-density 
are denoted j and p c , respectively. Note that the absence of a GW source terms in 
Eqs. (2.7) and (2.8) is valid only within the given approximation, as can be seen from 
e.g. Ref. [45]], where source terms that are fully nonlinear in the GW-amplitude 
are included. In addition to the explicit source terms in (2.9)-(2.10), naturally the 
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gravitational effects that is associated with the tetrad (2.3)-(2.4) must also be kept 
in mind. 

Next we turn our attention to the particle description. The equation of motion 
for a particle of mass m and charge q in an electromagnetic and gravitational wave 
field is 



-p = <Z [E+(7m)- 1 pxB] -G 



(2.11) 



where 7 = PiP i / (mc) 2 and the four-momenta is p a = jmdx a /dt. The gravi- 

tational force like term G % = F t ab p a p b /7m, where T l ab are the Ricci rotation coeffi- 
cients, becomes 



Gi = i(w ph - Pz/im) h+pi + h x p 2 
G 2 = ^(v ph -Pzhrn) -h+p 2 + hxPi 
G3 = |(7m) _1 h + (pl - pi) + 2h x p 1 p 2 



(2.12) 
(2.13) 
(2.14) 



for weak gravitational waves propagating in the z-direction in Minkowski space, 
where v p h is the phase velocity of the gravitational wave, which we at this point 
allow to deviate slightly from c. 

Next we apply kinetic plasma theory, representing each particle species by a 
distribution function / governed by the Vlasov equation. In tetrad form the Vlasov 
equation reads [59]] 

Cf = 

where the Liouville operator is 

C ee dt + (c/pV + [ f em - r\ bP a p h c/p a ] d pi 

and the electromagnetic force responsible for geodesic deviation is Fg M = q{E l + 
e % i k pjBk/ r )m). In vector notation the Vlasov equation reads 



d t f + 



P V/ 

7m 



q E+ 



p x B 

7m 



V P / = 



(2.15) 



where V p = (d Pl , d P2 ,d P3 ). In the absence of gravitational waves, the Vlasov equa- 
tion has the following spatially homogeneous (thermodynamical) equilibrium solu- 
tion, the Synge-Jiittner distribution, e.g. [40], 

n H __ u ^ 



fs.J 



(2.16) 



where no is the spatial particle number density, fi = mc 2 /kBT 7 ks is the Boltzmann 
constant, T the temperature and K 2 {n) is a modified Bcsscl function of second kind. 



3. Wave-wave interaction 

Next we let all quantities consist of a superposition of three waves of different kinds. 
Firstly we have a gravitational wave (frequency and wavevector (w g ,k g )), next an 
electromagnetic wave (w em , kem), and finally an electrostatic wave (Langmuir wave) 
(u> es , k cs ). Since we are dealing with high-frequency waves, the mass and charge will 
in what follows refer to electrons. The three waves are assumed to obey the following 
matching conditions 

LO g = LO cm + LO cs (3.1) 
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and 



kg — kgQs — k em + k cs (3-2) 

Since there is no external magnetic field within our model, the gravitational wave 
does not induce an electric field at the frequency and wavevector (w g , k g ). Further- 
more, the metric perturbations at the combinations (^1,2, ki.,2) can be neglected, 
as a consequence of the high-frequency approximation. 

As a prerequisite to the nonlinear calculations we first consider linear theory. Lin- 
earizing the Vlasov equation, in the absence of EM- fields, the perturbed distribution 
function of the GW is given by. 

_ 2ip a G a dfsj . . 

u s d\ P \ 2 

where we have introduced the notation u> g = uo g — k g cp 3 /p - Noting that the different 
GW-polarizations obey 

(uj 2 g - k 2 g c 2 ) h x = -2kST 12 = -2k J ^f e d 3 p (3.4) 

and 

(w* - k 2 g c 2 ) h+ = -2k(ST u - ST 22 ) = -2k J ^-ZM f gd 3 p (3.5) 
we deduce the same dispersion relation for both GW-polarizations, namely 

D g (u M = co 2 g - k\c 2 - 2ku J ^ffpP = (3.6) 

where we here and from now on use the normalization J fsjd 3 p = 1 of the un- 
perturbed distribution. In the high-frequency approximation, the last term of (3.6) 
is a small correction, comparable to the contribution from the background curva- 
ture [5]], and thus we may use the the approximation D g = u> 2 — k 2 c 2 . Similarly, 
from the linearized Vlasov equation (without gravitational fields), together with 
Maxwell's equations we deduce the dispersion relation for electromagnetic waves 

^ cm K m , fc cm ) = 1 - %! + ^- / ^^d 3 p = (3.7) 

Wem Wem J 7"cm 8\p\ 2 



2m ^P f k cs • p dfsj j3 



and for electrostatic waves 

k 2 s J JUJ CS d\p\ 2 

where u> p = (n q 2 j 'e^m) 1 / 2 is the electron plasma frequency. When nonlinear in- 
teractions are taken into account, the wave amplitudes will be time-dependent. We 
note that as far as the linear terms are concerned the only modification needed is 
the simple substitution D es E es = (dD es /div es )dE es /dt, where the tilde denotes the 
weakly time-dependent amplitude, and similarly for the other waves [57]]. Next, 
for definitcness, we assume the wave- vectors to span the plane perpendicular to e 2 . 
Furthermore, we note that for symmetry reasons, the h x -polarization in this geom- 
etry couples to the EM-wave polarized with the electric field along e 2 (nonlinearly 
combined with the electrostatic wave) , whereas the /i+ -polarization couples to the 
EM- waves with magnetic field along e 2 . These two cases are similar, and from now 
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on we limit ourselves to the former choice of polarization. The amplitude evolution 
for h x is found by keeping resonant terms proportional to both the EM and electro- 
static wave in the energy-momentum tensor. Next we use Eq. (3.4), generalized to 
keep the nonlinear terms and a weakly time-dependent amplitude. Substituting the 
expression for f g from the Vlasov equation, but now with second order nonlinear 
terms included, we obtain 



dh x 
s dt 

where the coupling coefficient C g is 

a 



2eonC„E cs E c 



(3.9) 



k cs • ei 


f V\V\ 




h k 
^es ^es 


J -fU) cm U) es 




1 k cm • k cs 

' \ 2 




dfsj 






d\p\ 2 



4p • k c 



d 2 fsj 
d(\p\ 2 ) 2 



d 3 p. 



(3.10) 



To find the EM-wave evolution, we now include all resonant nonlinear source terms 
for the EM-wave (effective gravitational currents, nonlinear terms involving the 
gravitational force, and nonlinearities coming directly from the tetrad), and solve 
for dE em /dt, in which case we obtain 



dD cm dE c „ 

dLUcm dt 



(3.11) 



where 



C e 



k 2 r 2 
"cm 1 - 



2uj c 



+ 



7.1 



vl 



7W om CJ C s 



■ei 2 



df. 



SJ 



d\ P [ 



4piw g d) em 



d(\ P \y 



d 3 p. (3.12) 



Finally the calculation is completed by doing analog calculation for the electrostatic 
wave evolution, and the result is 



where 





0D es dE cs _ 




du) es dt 


_ ^4 


'k cm • ei f p\ ( 


&cs 


2ui 2 J ju> es \ 


K - 


P3k s /^m e ) k es • ei + 


CJ cm CJ cs 



■CeshxE*^ 



(3.13) 



ei 



2k„ 



1 uj e pi 



ryuj cs uj g u; CI[L m e 



(2pi/7m e )k cs • k g ^ df S j 
d\p\ 



rd 3 p 



(3.14) 



The coupling coefficients C cs and C cm (although not C g ) has been calculated in Ref. 
[40]], using a coordinate frame formalism. As can be seen the previously computed 
coefficients arc algebraically much more complicated, and thus a comparison is 
difficult in general. For the special case where the temperature of fsj approaches 
zero, the coefficients simplify a lot, and we obtain 

C cs — C cm = — t ~ (3.15) 
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which agrees with Rcf. [40]]. Furthermore, we find that in our case that also the 
gravitational case coefficient is similar, and we can define the common coupling 
coefficient C = C cs — C om = C g . As a result of the agreement of the coefficients, 
the energy change of each of the waves in the cold limit can be written 

dW e 



dt 
dt 

and 

dW m 



= -u3gV, (3.16) 
= u cs V (3.17) 

= UcmV (3.18) 



dt 

where the energy density of each wave is 

W g =uj 2 s \h x \ 2 /2K, (3.19) 

Wc S = £ u es (dD es /duj es )\E es \ 2 , (3.20) 

Wc m (dD cm /duj cnl )\E cm \ 2 (3.21) 

and 

V = e Ch x E; m E; s + c.c (3.22) 

where c.c. denotes complex conjugate. The equations (3.16)-(3.18) together with 
(3.1) thus shows that the three- wave interaction process conserves the total wave- 
energy. 



4. Summary and conclusion 

In the present paper we have reconsidered the process of three-wave interaction 
between gravitational-, electromagnetic- and electrostatic waves in a collisionfrec 
plasma. Using an ortonormal tetrad description, the algebraic complexity of the 
coupling coefficients is much reduced, as compared to previous authors [40]]. Two 
of the coupling coefficients (for the electrostatic and the electromagnetic wave) can 
be shown to agree with previous results in the cold limit. The third coefficient (for 
the GW) has not been calculated before. Inclusion of the back-reaction on the GW 
makes it possible to deduce and energy conservation law (Eqs (3.16)-(3.18) together 
with (3.1)) in the cold limit. The energy conservation law is a natural consequence 
of the Manley-Rowe relations [57]], which makes the same coefficient appear in all 
three coupled equations (c.f. C = C cs = C om = C g .). However, this relation was 
only verified in the cold limit. Although the present coefficients were not extremely 
complicated, it has proved to be difficult to deduce whether the coupling coefficients 
are symmetric (i.e. fulfilling the Manley-Rowe relations), for the general case of a 
finite temperature. An interesting question, that remains for future research, is thus 
to decide whether this is due to some principal lack of a canonical Hamiltonian 
structure (that is the basic source of symmetric coupling coefficients, see e.g. [60]]) 
of the Einstein-Maxwell- Vlasov system, or if the difficulties are merely due to the 
algebraical complexity. 
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